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1.  INTRODUCTION 


The  description  of  charge  trapping  in  Si02  insulating  layers  in  a 
total  dose  environment  and  the  accompanying  device  threshold  voltage  shift 
has  been  of  long  standing  concern. The  accumulation  of  trapped 
positive  charge  in  gate  oxides  is  known  to  contribute  to  threshold  voltage 
shifts  of  a  few  tenths  of  volts  to  many  tens  of  volts  in  devices  and  test 
structures . ^ Because  the  threshold  voltage  shift  depends  on  the  first 
moment  of  the  oxide  charge  distribution,^^  the  spatial  dependence  of  this 
distribution  is  of  critical  importance.  A  number  of  experimentally 
accessible  factors  have  been  identified  which  contribute  to  the  shape  of 
the  oxide  charge  distribution.  Among  the  most  important  factors  are:^"^ 

(1)  total  dose  absorbed,  (2)  applied  field,  (3)  oxide  thickness,  and  (H) 
temperature.  The  physical  mechanisms  documented  as  relevant  for  a 
description  of  the  spatial  distribution  of  trapped  charge  in  a  radiation 
environment  include:  (1)  radiation  induced  electron-hole  pair  genera¬ 
tion, (2)  geminate  recombination,^"^' (3)  carrier  sweep 
out, (M)  hole  capture, (5)  electron  capture, ^ and 
tunneling. ^ 

In  the  following  section  we  develop  a  description  for  the  oxide- 
trapped  charge  distribution  in  a  low  dose  rate  environment.  The  spatial 
dependence  of  this  distribution  is  shown  to  be  a  consequence  of  the  spatial 
dependence  of  the  steady-state,  radiation-generated  carrier  distributions 
available  for  trapping.  The  steady-state  distribution  of  electrons  and 
holes  is  due  to  a  balance  between  depletion  of  these  mobile  carriers  in  the 
oxide  resulting  from  field  assisted  sweep  out  and  the  accumulation  of 
radiation-generated  electrons  and  holes  escaping  geminate  recombination. 

For  positive  gate  bias,  the  electron  flux  is  zero  at  the  oxide-semicon¬ 
ductor  interface  and  increases  toward  the  gate.  The  hole  flux  is  zero  at 
the  gate  and  increases  toward  the  ox ide-semi conductor  interface.  Hie 
resulting  trapped  charge  distribution,  due  to  electron  and  hole  trapping, 
is  shown  to  have  a  maximum  near  the  oxide-semiconductor  interface  for 
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positive  bias,  even  for  a  constant  trap  distribution.  The  trapped  charge 
distribution  is  integrated  to  yield  an  expression  for  the  trapped  oxide 
charge  contribution  to  the  threshold  voltage  shift.  Two  experimentally 
important  quantities  are  derived:  (1)  the  initial  midgap  voltage  shift  per 
dose  and  (2)  the  maximum  (steady  state)  midgap  voltage  shift.  The  field 
dependence  of  these  quantities  is  shown  to  be  a  consequence  of  the  field 
dependence  of  the  geminate  recombination  escape  probability  and  the  field 
dependence  of  the  electron  and  hole  capture  cross  sections.  The  change  in 
the  trapped  charge  distribution  due  to  tunneling  at  the  oxide-semiconductor 
interface  is  derived  and  compared  to  the  results  without  tunneling.  The 
effect  of  internal  fields  as  a  consequence  of  charge  trapping  is  derived 
and  compared  to  the  results  without  correction  for  the  internal  field. 

In  the  subsequent  section,  experimental  results  for  the  initial  midgap 
voltage  shift  per  dose  and  the  maximum  midgap  voltage  as  a  function  of 
field  are  compared  to  the  theoretical  results. 
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1 1 .  THEORY 


The  rate  equation  for  hole  trapping  in  Si02  niay  be  developed  by 
considering  the  continuity  equations  for  radiation-generated  valence  band 
holes,  conduction  band  electrons,  and  the  continuity  equation  for  trapped 
holes  separately.  Electron  trapping  is  assumed  to  take  place  by  a  positive 
coulombic  trapping  mechanism  and  hole  trapping  by  neutral  traps.  The  one¬ 
dimensional  continuity  equation  for  conduction  electrons  is 


3n  ^--^n  • 

3t  ’  3x  "o^*^  ■  “n'^n^T 


where 


the  conduction  electron  concentration 
electron  flux 

electron-hole  pairs  generated  per  dose 
dose  rate 

electron  capture  cross  section 
(coulombic  traps) 

trapped-hole  concentration 

electron-hole  escape  probability. 


Similarily,  the  continuity  equation  for  the  valence  holes  is 


S  =  ^  -  '>pJp<“T  -  '“t> 
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where 


P  = 

°P  ^ 

Jp  = 

The  continuity  equation  for  trapped  holes  is 

gp 

IT  =  «pJp'“T  -  fT*  -  ”nJn'’T-  <3> 

For  typical  values  of  capture  cross  sections  and  trapped  charge 
concentrations  the  radiation  generation  term  in  Eqs.  (1)  and  (2)  domi¬ 
nate.  Under  positive  gate  bias,  holes  are  swept  toward  the  interface, 
making  the  concentration  of  holes  available  for  trapping  zero  at  the  gate, 
while  electrons  are  swept  toward  the  gate,  making  the  concentration  of 
electrons  zero  at  the  oxide-semiconductor  interface.  In  a  low  dose  rate, 
high  total  dose  environment,  we  assume  that  the  conduction  electrons  and 
valence  band  holes  reach  steady  state  in  a  short  time  compared  to  the  time 
for  irradiation.  For  electrons  this  occurs  within  a  fraction  of  a  pico¬ 
second  for  oxides  less  than  1000  A.  For  holes  steady  state  will  be  reached 
in  a  fraction  of  a  microsecond.  As  our  irradiations  lasted  from  tens  of 
minutes  to  a  few  hours,  the  assumption  of  a  steady  state  distribution  is 
reasonable.  Under  these  conditions  Eqs.  (1)  and  (2)  may  be  solved  to  yield 

Jp  ' 


concentration  of  oxide  hole  traps 
valence  band  hole  concentration 
hole  capture  cross  section  (neutral  traps) 
hole  flux. 


Jp  =  n^D^x 


(ilb) 
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where  is  the  oxide  thickness.  The  origin  of  the  coordinates  is  at  the 
gate  electrode.  The  rate  equation  for  hole  trapping  is  then 


dP 

dT  =  -  ’“t*  -  °n"o^*<‘ox 


X)P,p. 


This  equation  may  be  solved  to  yield: 


=  N,pf(x)(1 


g-T(x)"^t 


) 


where 


f(x)  =  o  x/(a  X  +  o  (t  -  x)) 
p  p  n  ox 


(5) 


•  (6a) 


(6b) 


t(x)~^  =  n  D()!(o  X  +  0  (t  -  x)),  (6c) 

op  n  ox 

From  this  point  on  we  assume  that  is  constant  unless  specifically  stated 

otherwise.  Note,  Eq.  (6a)  predicts  a  spatially  dependent  trapped  charge 

distribution,  even  for  constant  N.p.  The  spatial  variation  of  the  trapped 

charge  distribution  results  from  the  spatial  variation  of  the  available 

carriers  as  given  by  Eqs.  (4a)  and  (4b)  and  the  electron-’  'e  trapping 

kinetics  as  described  by  the  electron  and  hole  capti'ie  cross  sectx.  , 

and  0  . 

P 

To  describe  the  effect  of  this  distribution  on  the  midgap  voltage,  the 
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first  moment  of  the  charge  distribution  must  be  evaluated,  thus  ' 


qN_  ^ox 

AV  =  —  r  X  f(x)(1 
^  0 


(7) 


where  AV  is  the  magnitude  of  the  midgap  voltage  shift,  q  is  the  electronic 
charge,  and  e  the  oxide  permittivity. 
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Two  quantities  of  experimental  interest  are  the  initial  rate  of  change 
of  the  midgap  voltage  and  the  maximum  midgap  shift.  Upon  substitution  of 
Eq.  (6)  and  expansion  for  small  t,  integration  of  Eq.  (7)  yields  the 
initial  midgap  voltage  shift  per  dose,  thus 


AV/D 


t-H3 


i—  N_t^  n  0  0 
3e  T  ox  o  p 


(8) 


where  D  =  Dt .  This  quantity  is  directly  proportional  to  the  number  of 
available  traps,  the  generation  rate  of  electron-hole  pairs,  the  fraction 
of  electrons  and  holes  which  escape  recombination,  and  the  hole  capture 
cross  section.  Electron  trapping  does  not  affect  the  initial  trapping  rate 
because  there  are  too  few  trapped  holes  to  accommodate  any  electrons. 


The  maximum  voltage  shift  is  obtained  by  setting  the  exponential  term 
in  Eq.  (7)  to  zero.  Substitution  of  Eq.  (6b)  and  integration  yields 


AV 


max 


(9a) 


where 


N_(o  /o  )  .  P 

^ - _  {£n(o  /o  )  +  o(On^°n  ■ 

(a  /o  -  1)3  P  "  2  p  n 

P  n 


eff 


(Oo/c^  -  1)}. 


(9b) 


The  maximum  midgap  voltage  shift  is  critically  dependent  on  the  hole-to- 
electron  capture  cross  section  ratio.  In  the  limit  that  the  electron 
capture  cross  section  goes  to  zero  (no  electron  trapping)  reduces  to 

Wj./2  and  the  effect  of  either  capture  cross  section  is  lost. 

As  seen  in  Eq.  (8)  any  electric  field  dependence  of  the  initial  slope 
must  be  attributed  to  the  field  dependence  of  the  hole  capture  cross 
section,  Op,  and  escape  probability,  ♦,  all  other  quantities  being  field 
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independent.  The  field  dependence  of  the  maximum  midgap  voltage  shift  is 
due  to  the  field  dependence  of  the  ratio  of  the  hole  to  electron  capture 
cross  sections.  The  field  dependence  of  the  capture  cross  sections  and  the 
escape  probability  are  discussed  in  the  next  section, 

A.  FIELD  DEPENDENCE 

We  assume  that  the  escape  probability  is  dominated  by  geminate 
recombination.^®  For  small  field,  the  escape  probability  is  given  by^^ 

it(E)  =  4,^(1  +  qE/8iTe(kT/q)^)  (10) 

where  k  is  the  Boltzmann  constant,  T  is  the  absolute  temperature,  and  4^  is 
the  zero  field  escape  probability.  The  escape  probability  increases 
linearly  with  field  from  4^  to  unity.  For  larger  fields  the  escape 
probability  remains  unity  as  all  of  the  electron>hole  pairs  produced  by 
radiation  are  separated.  Shown  in  Fig.  1  is  the  escape  probability,  41(E), 
versus  the  voltage  across  a  225  A  oxide  at  room  temperature,  assuming  the 
escape  probability  is  unity  at  about  one  volt. 

We  may  describe  the  field  dependence  of  the  electron  capture  cross 

20  2 1 

section  by  appeal  to  a  cascade  capture  model  of  a  coulombic  trap.  »  In 
zero  applied  field  the  potential  of  the  trap  is 

V(r )  =  -q/4ii£r.  (11) 

The  approaching  electron  is  captured  when  this  potential  is  equal  to 
-2kT/q,  That  the  maximum  probability  of  capture  occurs  at  a  potential  of 
-2kT/q  has  been  shown  to  be  a  consequence  of  the  cascade  capture  of 
carriers. If  the  trap  is  deep  enough  the  energy  levels  near  the  top  of 
the  trap  are  closely  spaced.  The  captured  carrier  is  initially  trapped  in 
an  energy  level  consistent  with  the  energy  dissipation  mechanism  in  the 
material,  the  most  probable  of  which  is  single-phonon  creation  at  energies 
comparable  to  kT.  After  the  ir.itial  capture  the  trapped  carrier  cascades 
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to  lower  ene-^gy  levels.  Thus,  the  zero  field  electron  capture  cross 
section  may  be  estimated  as 


o 


no 


where 


r^  =  q/8iie(kT/q) . 


(12a) 


(12b) 


Fig.  1.  Geminate  Recombination  Escape  Probability, 

Reduced  Electron  Capture  Cross  Section, 
and  Reduced  Hole  Capture  Cross  Section  vs 
Potential  Across  a  225  A  Oxide 

At  room  temperature  this  gives  a  =  1.6  *  10"^^  cm^  which  is  consistent 

no 

with  previously  cited  values.®’ In  an  applied  field  the  electron 
capture  cross  section  may  be  calculated  assuming  Frenkel-Poole  barrier 
lowering. The  perturbed  potential  in  the  direction  of  the  applied  field 
is : 
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(13) 


V(x)  =  -(1/4iiex)  -Ex 

o 

where  is  the  magnitude  of  the  macroscopic  field.  To  solve  for  the  field 
dependent  capture  radius,  r',  we  require  that  the  potential,  \/(r),  at  r'  be 
2kT/q  lower  than  the  extremum  of  V(x).  The  result  is 

r'  =  (x'  +  x'^/2r^)  -  ((x'  ^  x'^/2r^)  -  (1i<a) 

where  x'  is  the  value  of  x  in  Eq.  (13)  at  the  extremum: 

X'  =  (14b) 

In  Fig.  1,  a  /a  is  shown  as  a  function  of  voltage  across  a  225  A  oxide  at 
n  no  ° 

room  temperature. 

The  field  dependence  of  the  hole  capture  cross  section  is  calculated 
in  much  the  same  way,  except  that  the  hole  trap  is  neutral.  In  zero 
applied  field  the  neutral  trap  is  polarized  by  the  field  of  an  approaching 
hole.  The  induced  dipole  potential  of  the  trap  is 

V(r)  =  -qR^/47ter^  (15) 

where  Rq  is  the  average  radius  of  the  trap  (atomic  dimensions).  For  R^ 

about  equal  to  3  A,  o  =  1.4  x  10“^^  cm^  at  room  temperature  and  is 

1  7  fi  in  ?4 

consistent  with  previously  cited  hole  capture  cross  sections.  ^ 

Frenkel-Poole  barrier  lowering  and  cascade  capture  leads  to  a  fifth  order 
equation  for  the  field  dependent  capture  radius.  Over  the  applied  fields 
of  interest  here  (less  than  4  x  10^  V/cm)  the  field  dependent  hole  capture 
cross  section  is  closely  approximated  by 

0  =  o_(1.0  +  1.9  X  10"^  e:^^)"’  (16) 

p  po  o 
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where  has  units  of  V/cm.  °p^°pQ  ®  function  of  applied  voltage  across 

a  225  A  oxide  at  room  temperature  is  shown  in  Fig.  1,  Comparison  of 

and  a^/a  shows  that  r  ''r  is  much  more  sensitive  to  the  applied  field, 
n  no  n  no  rr  » 

particularly  at  lower  fields.  Within  .5  volts  (0  -  2  x  10^  V/cm) 

drops  by  more  than  an  order  of  magnitude  while  o  /Op^  drops  by  less  than 

^5%. 


B.  TUNNELING 


Tunneling  of  holes  out  of  the  oxide  (or  equivalently,  electrons  into 
the  oxide)  at  the  oxide-semiconductor  interface  is  known  to  be  an  important 
phenomenon  in  determining  threshold  voltage  shifts  in  MOS  devices. ^ 

We  have  included  tunneling  in  our  development  by  adding  a  tunneling  term  to 
Eq.  (3)  of  the  form 


Tunneling 


-P.j,g(x,E,t) 


(17) 


where  g(x,E,t)  is  the  tunneling  rate.  For  small  t,  Pj  s  0,  tunneling  is 
unimportant  and  the  results  for  the  initial  rate  of  change  of  the  midgap 
voltage  [Eq.  (8j  remain  unchanged.  At  saturation  the  maximum  midgap  volt¬ 
age  shift  will  change  due  to  tunneling  because  the  trapped  hole  distribu¬ 
tion,  f(x),  will  change.  The  new  form  of  the  trapped  charge  distribution 
which  includes  tunneling  is 


f.p(x)  =  OpX/(OpX  +  -  x)  +  g(x,E)/n^D(t.) .  (I8) 

In  the  simplest  form  of  the  WKB  approximation,  the  tunneling  rate  g(x,E)  is 

g(x,E)  =  g^e^^^^^ox”*^  (19) 


where 


8  = 


(19b) 
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m  is  the  tunneling  effective  mass,  the  trap  depth,  h  is  Planck's 
constant  divided  by  2n ,  and  is  the  attempt  frequency.  The  attempt 
frequency  may  be  estimated  with  the  aid  of  the  Heisenberg  uncertainty 
relations 

g  >  h/l6m«o  .  (20) 

o  p 

Shown  in  Fig.  2  is  a  comparison  of  the  saturation  charge  distributions  with 
and  without  tunneling  versus  x/t^jj  for  two  applied  fields  (gate  potential 
positive).  The  tunneling  probability  was  evaluated  with  m*  equal  to 
.42  free  electron  masses  and  equal  to  3-1  eV.  Without  tunneling  the 
trapped  charge  distribution  function,  f(x),  rises  steadily  from  zero  at  the 


Reduced  Oxide  Position 


gate  to  unity  at  the  semiconductor  interface.  With  tunneling  the  charge 
distribution  function,  fj(x),  drops  abruptly  near  the  semiconductor  inter¬ 
face.  In  a  very  short  distance  (less  than  10  A  for  a  225  A  oxide)  the 
charge  density  decreases  by  two  orders  of  magnitude.  Because  the  maximum 
midgap  voltage  shift  depends  on  the  first  moment  of  the  charge  distribution 
the  loss  of  charge  near  the  semiconductor  interface  due  to  tunneling  will 
have  a  large  effect  on 

The  Ngff  which  accounts  for  tunneling  is: 


eff 


N_(o  /o  ) 

- — ^ - 5  {an((o  /o  -  1)(1  -  Ax/t  )  +  1) 

3  ‘  p  n  ox 


(o  /o  -  1 ) 
P  n 


.  -  1)^  (1  - 


(V°n  -  ')<'  -  ‘’•'‘ox” 


(21) 


where  Ax,  the  thickness  of  the  charge  depleted  layer  due  to  tunneling  is 


Ax  s  (2e)*^tn[g^/n^D4.t^^(Op  +  0.2o^)].  (22) 


C.  INTERNAL  FIELDS 

Due  to  the  distribution  of  trapped  charge  shown  in  Fig.  2  the  field 
internal  to  the  oxide  will  be  different  from  the  applied  field.  »  ’ 
Using  Poisson's  equation  we  have  calculated  the  average  field  in  the 
trapped  charge  region  of  the  oxide.  With  the  boundary  condition  of 
positive  applied  gate  bias  during  irradiation,  the  result  of  this 
calculation  is: 


qN_.Ax 


E  =  E  -  - 7— 
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(23) 


Superimposed  on  Fig.  2  is  the  charge  distribution  fj(x)  evaluated  at  E.  The 
result  of  the  reduced  internal  field  at  the  applied  fields  shown  reduces 
the  trapped  charge  distribution. 

In  the  next  section  the  analytic  results  derived  above  are  compared  to 
experiment . 


15 


III.  EXPERIMENTAL  COMPARISON 


Test  capacitors  used  in  this  study  consisted  of  a  thermal  225  A  oxide 
grown  on  p-type  silicon,  with  a  tungsten  silicide/n+  polysilicon  gate.  The 
capacitors  were  annealed  at  lOOCC  in  nitrogen  for  15  minutes.  The  test 
capacitors  were  irradiated  under  bias  at  room  temperature  using  a  Shepherd 
Model  81-22  Co-60.  All  dosimetry  and  sample  irradiations  were  done  in  a 
lead/aluminum  box  to  limit  dose  enhancement  effects  due  to  cavity  geometry. 
The  radiation  response  of  the  MOS  capacitors  was  monitored  in  situ  using 
high  frequency  (1  MHz)  capacitance-voltage  measurements.  The  dose  rate  in 
all  cases  was  156  kRad(Si)/hour.  The  threshold  voltage  shift  due  to 
trapped  oxide  charge  was  determined  as  the  midgap  shift  of  the  post- 
irradiation  C-V  measurements.  Interface  state  generation  was  negligible  as 
determined  by  C-V  stretchout  measurements. 

Shown  in  Fig.  3  is  a  typical  result  for  the  midgap  shift  versus  total 
dose  for  a  given  applied  field.  In  Figs.  3-5,  AV  denotes  the  magnitude  of 
the  midap  voltage  shift  and  therefore  is  plotted  along  the  positive  axis 
even  though  the  measured  shifts  were  negative.  The  uncertainty  in  the 
voltage  shift  is  ±10  mV.  The  solid  curve  is  an  empirical  fit  to  the  data 
to  aid  in  determining  the  initial  slope.  Nine  such  sets  of  data  were 
obtained  for  applied  oxide  fields  from  zero  to  4.2  x  10^  V/cm  and  the 
initial  voltage  shift  per  dose  (initial  slope)  determined.  These  results 
are  shown  in  Fig.  4.  The  heavy  dots  are  the  values  for  the  initial  voltage 
shifts  per  dose  as  determined  by  the  empirical  fit.  For  comparison,  the 
average  shifts  per  dose  at  50  krad(Si}  are  shown  as  crosses.  The  error 
bars  represent  the  uncertainty  in  determining  the  initial  slope.  The 
horizontal  scale  is  given  as  voltage  across  the  oxide. 


TOTAL  DOSE  (rad  Si))  10^® 

Fig.  3.  Typical  Data:  Magnitude  of  the  Midgap 

Voltage  Shift  vs  Dose  Data  (Vg  =  '  volt). 

Superimposed  on  Fig.  4  are  the  results  obtained  by  using  Eqs.  (8)  and 
(16)  of  the  previous  section.  We  assumed  a  value  of  4  x  10®  cm“^  for 
UpOp^,  which  is  very  close  to  previously  cited  values  for  this  quantity,^^ 
and  a  zero  field  escape  probability,  of  1.5%  which  is  also  very  close 
to  previously  cited  values.^  As  seen  the  data  are  well  presented  by  the 
theory.  The  low  field  rise  in  the  initial  midgap  voltage  shift  per  dose  is 
dominated  by  the  linear  increase  with  field  of  the  geminate  recombination 
escape  probability.  At  high  fields  (above  4.4  x  10®  V/cm)  the  escape 
probability  saturates  at  unity  and  the  field  dependence  is  dominated  by  the 
decrease  of  the  hole  capture  cross  section  with  field.  For  a  more  real¬ 
istic  form  of  the  geminate  recombination  escape  probability  (i.e.,  a  less 
abrupt  transition  to  saturation)  the  cusp  shown  in  Fig.  4  would  give  way  to 
a  smoother  low  to  high  field  transition.  Examination  of  Eq.  (16)  shows 

that  the  high  field  behavior  of  the  midgap  shift  per  dose  should  decrease 
-  55 

as  •  This  high  field  behavior  has  been  seen  experimentally 
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elsewhere.  We  are  unaware  of  any  theoretical  Justification  for  this 
behavior  prior  to  our  analysis  of  the  previous  section. 


Fig.  4.  Magnitude  of  the  Initial  Midgap  Voltage  Shift 
per  Dose  vs  Potential  Across  a  225  A  Oxide 


Shown  in  Fig.  5  is  versus  the  applied  voltage,  V^,  across  the 

oxide.  The  dashed  curve  is  the  theoretical  result,  Eq.  (9),  with  no 
tunneling  and  no  correction  due  to  the  internal  field.  The  dash-dot  curve 
is  the  result  for  in  which  tunneling  is  included  using  Eq.  (21)  for 

Neff  and  assuming  that  Ax,  the  thickness  of  the  depleted  layer,  is  50  A. 
The  result  represented  by  the  dash-dot  curve  excludes  the  internal  field 
correction.  As  seen,  tunneling  dramatically  decreases  the  predicted  maxi¬ 
mum  midgap  voltage  shift  because  charge  is  removed  from  the  peak  of  the 
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distribution  which  is  close  to  the  oxide-gemiconductor  interface.  The 
solid  curve  is  the  theoretical  result  predicted  by  including  tunneling  and 
the  correction  of  the  internal  field  Eq.  (23).  Because  the  average 
internal  field  over  the  oxide  charge  distribution  is  smaller  than  the 
applied  field,  less  charge  is  trapped  (see  Fig.  2).  Thus,  the  maximum 
midgap  voltage  shift  is  reduced.  Much  below  2  volts  (-  9  *  10^  V/cm)  our 
fomulation  of  the  Internal  field  correction  leads  to  a  negative  average 
field  in  the  oxide,  which  violates  the  original  assumptions.  Therefore,  we 
neglect  to  show  any  theoretical  results  below  about  2  volts.  The  reason 
for  this  breakdown  is  that  the  charge  distribution  function,  f<p(x),  in  the 
correction  term,  Eq.  (23) >  has  been  evaluated  at  the  applied  field.  At  low 
fields  the  difference  between  the  charge  distribution  evaluated  at  the 
applied  field  and  the  internal  field  is  relatively  large  (see  Fig.  2). 
Therefore,  the  correction  is  too  large,  leading  to  a  breakdown  of  the  first 
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order  model  at  low  fields.  At  high  fields,  the  difference  in  the  distribu¬ 
tion  function  at  the  applied  field  compared  to  the  internal  field  is  slight 
(see  Fig.  2).  Higher  order  corrections  to  the  internal  field  model  will, 
therefore,  preserve  the  high  field  results  and  improve  the  low  field 
results. 

Superimposed  on  Fig.  5  are  the  maximum  midgap  voltage  shifts  as 

determined  from  our  data.  The  triangles  represent  AV  determined  from 

max 

our  measurements  by  extrapolation  to  infinite  dose.  All  measurements  were 
terminated  at  1  *  10^  rad  (Si).  The  error  bars  reflect  the  uncertainty  in 
AVmax  to  the  extrapolation.  In  the  region  above  2  volts  the  theory 
predicts  the  data  reasonably  well.  Due  to  tunneling,  the  effect  of  the 
internal  field,  and  the  balancing  of  electron  and  hole  trapping,  the 
maximum  midgap  voltage  shift  as  the  applied  field  is  increased  is  greatly 
reduced  compared  to  the  results  predicted  by  ignoring  these  effects. 
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IV.  DISCUSSION 


We  have  derived  expressions  for  the  initial  slope  of  the  midgap 
voltage  shift  and  the  maximum  midgap  voltage  shift  as  functions  of  the 
applied  field  across  the  oxide  in  a  low  dose  rate  environment.  Electron 
trapping  as  well  as  hole  trapping  is  included  in  the  formulation.  Hole 
trapping  in  the  oxide  takes  place  by  cascade  capture  by  polarized  neutral 
traps.  Electron  capture  takes  place  via  coulombic  cascade  capture  by 
trapped  holes.  The  distribution  of  holes  and  electrons  available  for 
capture  is  shown  to  be  a  steady  state  distribution  in  which  the  fraction  of 
holes  and  electrons  captured  is  small  compared  to  the  number  of  radiation 
generated  electron-hole  pairs  that  escape  geminate  recombination.  It  has 
been  shown  that  the  field  dependence  of  the  initial  slope  of  the  midgap 
voltage  shift  is  due  to  the  field  dependence  of  the  hole  capture  cross 
section  and  geminate  recombination  escape  probability.  The  maximum  midgap 
shift  is  shown  to  depend  on  the  ratio  of  the  hole  to  electron  capture  cross 
sections.  The  field  dependence  of  the  capture  cross  sections  were  derived 
using  Frenkel-Poole  barrier  lowering.  Tunneling  near  the  semiconductor 
interface  has  also  been  included.  We  have  also  included  an  Internal 
electric  field  due  to  the  trapped  charge  distribution  in  the  oxide. 

The  formulation  presented  provides  an  understanding  of  the  E 
decrease  in  the  high  field  behavior  of  the  initial  midgap  shift  per  dose 
cited  in  the  literature. The  theory  has  been  validated  by  comparison 
with  data  obtained  on  225  A  thermal  oxide  on  p-type  silicon  test  capacitors 
irradiated  under  bias  at  room  temperature.  The  theoretical  description 
explicitly  addresses  the  problem  of  determining  the  oxide  trapped  charge 
distribution  and  the  resulting  midgap  voltage  shift  as  a  function  of  total 
dose  and  applied  field. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  li.vestlgations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  Including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -field-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  LJttery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory;  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  mlcrowave/mlllimeter  wave  measurements, 
diagnostics  and  radioaetry,  aicrowave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  anu  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatu.cs  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospheric  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


